load('qingdaohyper.mat')
X1=double(qingdao11_03)';
X2=double(qingdao11_05)';
X3=double(qingdao11_06)';
X1=[X1(1:49,:);X1(51:93,:);X1(101:145,:);X1(148:149,:);X1(153:197,:)];
X2=[X2(1:49,:);X2(51:93,:);X2(101:145,:);X2(148:149,:);X2(153:197,:)];
X3=[X3(1:49,:);X3(51:93,:);X3(101:145,:);X3(148:149,:);X3(153:197,:)];
X1=Normalization(X1',1);
X2=Normalization(X2',1);
X3=Normalization(X3',1);
label=qingdao14_truthmap03;
label2=qingdao14_truthmap05;
label3=qingdao14_truthmap06;
clear qingdao* 
clear water tree stoneandsand  soil class* data_11 grass industry building
[r,c]=Corr_inf(X1,X2,label3,label2,par,graph,loc03,loc05);
  X1train=X1(c,:);
        X2train=X2(r,:);
    [Q, k, X1_mean, X2_mean] = Procrustes (X1train, X2train);